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Abstract 

Theories with extra dimensions have gained much interest in recent years as candidates for a 
possible extension of the SM. The observation of large extra dimensions through real graviton 
emission is one of the most popular related new phenomena. The main experimental signa- 
tures from graviton emission are production of single jet and single photon events, which have 
been studied in great detail. This work describes the implementation of graviton production 
together with either a Z or a photon in Pythia 8. The potential of using Z plus graviton 
production at the LHC as a complementary channel is also studied. For completeness, this 
work also includes the more recently proposed scenario of unparticle emission, since the ef- 
fective theory of unparticles to some extent represents a generalization of the large extra 
dimension model. 



1 Introduction 



The possibility of observing gravity in extra dimensions (ED) at the TeV scale was proposed a 
few years ago and has since gained a large interest. One of the first proposals was the so-called 
ADD scenario [1] which suggests that gravity alone would propagate in large extra dimensions 
(LED). This would explain the weakness of gravity experienced by the SM fields, since they 
only would have access to the normal 4 dimensional space-time. This scenario has been studied 
in detail [2, 3] and dedicated searches have been performed at different experiments [4, 5, 6]. 
The large extra dimensions in the ADD scenario have a flat space topology and, for this reason, 
astrophysical observations imply stringent constraints in the case of a small number of extra 
dimensions. It has, however, later been shown that these constraints would be evaded if space- 
time is slightly warped, as in the popular RS models [7]. In this case the graviton (G) would 
acquire a small effective mass, which makes the model [8] insensitive to astrophysical constraints 
from low-energy processes 1 . Then even one large extra dimension could have escaped the present 
experimental measurements [9]. 

Recently a so-called unparticle (U) model has also gained much attention [10, 11]. This 
relates to phenomena from a scale invariant sector which is coupled to the SM by a connector 
sector with a high mass scale. This scenario is normally considered to be less well motivated 
than extra dimensional gravity, which could solve the so-called hierarchy problem and is also 
motivated by string theory. It could, however, imply unusual experimental signatures at the LHC 
which should not be missed. In addition, the unparticle model is, from a phenomenological point 
of view, a generalization of the large extra dimension case and therefore both cases can to some 
extent be covered at the same time. 

The main experimental signatures from both gravity in large extra dimensions as well as 
the unparticle model would be an excess of single jet and single photon events. In this work we 
describe the implementation of Z/j + G/U production as a so-called semi-internal process to 
Pythia 8 [12]. Furthermore Z + G/U production at the LHC is studied, since this process could 
provide a complement in order to constrain the model if a signal would be observed in the main 
channels. 

The next section gives a reminder of the relevant model parameters and, in order to simplify 
comparison between different papers, the most common conventions are summarized. This is 
followed by the implementation of the process and the validation tests that have been performed. 
Finally a study at generator level of the pp — > Z + G/U cross section predicted at the LHC is 
presented, including a relative comparison with the similar SM process pp — > ZZ — > llvv. 

2 The Models and Parameters 

The extra dimension model considered in this work is the same as that described in [2]. Here all 
SM fields are confined to a 4 dimensional brane present in a larger dimensional space where only 
gravity can propagate. It is assumed that the brane is rigid so that effects from fields related 
to the brane dynamics can be neglected. The extra dimensions would be compactified with a 
radius R. In the case of real graviton production at the LHC the momentum component of the 
graviton in the extra dimensions would be observed on the brane as a mass. The finite size of 
the ED implies a discrete series of allowed mass modes, i.e. the so-called Kaluza-Klein (KK) 

lr Qie high-energy phenomena relevant at collider experiments are, however, unchanged. 
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tower. The graviton is coupled to the SM fields through the energy-momentum tensor as, 

AC = -(8irG N )^G$T^ (1) 

where Gn is the Newton constant and i is the KK mode index. For the graviton production 
cross section the extremely weak coupling to gravity is compensated by the additional phase 
space in the extra dimensions. For this reason gravitational phenomena would appear at energies 
around the D-dimensional fundamental scale of gravity, My, . The size of the ED in this scenario 
is allowed by experimental data to be almost as large as a millimeter [5]. Since the KK mode 
separation relates to the size as Am ~ 1/-R, the discrete KK series can be approximated by a 
continuous spectrum. In order to simplify the calculations, the extra dimensions are assumed to 
have the geometry of a n-dimensional torus. In this model of gravity in large extra dimensions, 
the size of the extra dimensions (R) and the fundamental scale of gravity (My>) are related to 
the Newton constant for n extra dimensions as, 

G" 1 = C ■ R n M™ +2 (2) 

where the exact definitions vary slightly between different papers. This work follows the pa- 
rameter definitions used in [2] where C = 8ir. However, to simplify for comparison between 
papers the different conventions are summarized. In the papers [2, 3, 13] the fundamental scale 
of gravity is referred to as My>,M and Ms- These are the three most common definitions and 
are related as follows, 

M n+2 = 2M£+ 2 (3) 
M™+ 2 = Stt 1 -? r (I) M£ +2 (4) 

The work in [13] also defines the size R slightly different compared to [2, 3], Rs = InRy. From 
these relations the different expressions in the literature for the integral over the very dense KK 
states reduce to, 

/_n/2 nil /™2\n/2-l 
dm2 ■ *„■(„,) (8) 

Here k is the graviton momentum component in the extra dimensions and m 2 is the graviton 
mass squared. 

The unparticle model studied here is described in [11]. In this scenario the so-called unpar- 
ticles originate from a scale invariant sector with a non-trivial fixed point. This sector interacts 
with the SM fields through a connector sector with a high mass scale, My. Renormalization 
effects in the scale invariant sector give rise to dimensional transmutation of the unparticles at 
an energy scale Kjj. This transmutation is determined by a scale dimension parameter dy which 
has to be greater than one due to unitarity arguments, but is allowed to take non-integer val- 
ues. This scenario gives rise to quite unusual phenomena. In terms of real unparticle emission 
it implies that the unparticle would not have a fixed invariant mass squared, but instead an 
invariant mass spectrum. The unparticle would also appear with a d[/-body final state phase 
space, which would give rise to an unusual missing energy signature. In accordance with [11] 
only the production of unparticles from the allowed spin-1 and spin-2 effective unparticle oper- 
ators is considered here and the unparticles are assumed to be SM singlets. In this model the 
cross sections are determined by the three unparticle parameters, dy, h-u and A. Here A is an 
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effective coupling between the unparticle and SM operators which is related to the connector 
scale M\j. In the case of spin-2 unparticles two different effective operators are allowed which 
are associated with the independent coupling constants A and A'. 

As pointed out in [11] the effective theory of spin-2 unparticle emission with A = A' is, 
apart from constant factors, equivalent to graviton emission in large extra dimensions. Here the 
invariant mass spectrum from the dimensional transmutation of the unparticles is identical to 
the dense KK tower of the large extra dimensions. For this reason the graviton emission process 
can be seen as a special case of the unparticle model and, as will be described below, both 
scenarios can be covered by the same implementation with only minor changes. This has the 
advantage of allowing the more motivated extra dimension search to be extended to cover the 
unparticle case. To cover both scenarios at once could be of great value in the case where future 
developments discover new and perhaps even more interesting analogies between the unparticle 
picture and gravity in extra dimensions. A potential connection between a non-integer dy value 
and so-called warped extra dimensions has for example been addressed in [14]. 



3 Process Implementation 

The process ff—>Z + G/U was implemented 2 as a semi-internal process to Pythia 8.108. It 
contains a 2-to-2 parton level process class which can be used by the main Pythia library. The 
class is implemented with the same structure as used for the internal parton level processes and 
is, therefore, also used in the same way as the internal classes inside the Pythia library. No 
Z/j* interference effects were taken into account and the Z decays isotropically. The location 
and details of the software can be found in the appendix. 
The differential cross section used corresponds to 

da (ff 7 .urn l>f| 2 A(du) fp 2 uY U ~\r 0w2. 

^ ff ^ z + u) = m^-2^7q;{7q;) e ^n Pu ) (6) 

where pfj is the invariant mass squared. Apart from the matrix element, the first term only 
contains standard 2-to-2 scattering phase space factors. The remaining terms contain the phase 
space and final state mass spectrum of the unparticle. Here, A{du) is a normalization constant of 
the dry-body phase space. The full spin- and color-averaged matrix elements given by equations 
(42) and (47) in [11] are used. This corresponds to the most general form of the matrix element 
from the two allowed effective spin-2 unparticle operators. The spin-1 case only includes the 
contribution from the vectorial operator. This is sufficient when considering unpolarized particle 
beams, since including the second axial-vectorial operator would result in the same formula but 
with A -► y/X* + X£. 

The variable mass spectrum shown in equation (6) gave rise to the main difference with 
respect to the internal processes available in Pythia. It turned out, however, to be conveniently 
implemented by re-weighting an internally produced Breit-Wigner distribution. By using the 
graviton already available in Pythia, the unparticle masses are generated according to a Breit- 
Wigner which is re-weighted to a flat distribution. The invariant mass measure of the cross 
section, {pjj) du ~ 2 dpjj, then ensures that the final events are generated with the correct mass 



2 As a by-product also the process // — » 7 + G/U was implemented in a separate class. This corresponds to 
the photon limit of the Z process as described later in the text. 
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distribution. In order to achieve a high MC efficiency one has to take care to match the bulk of 
the Breit-Wigner to the bulk of the unparticle distribution. It should be pointed out that this 
only affects the speed of generating events, by throwing events more often in the region where 
the differential cross section is large, and not the final cross section or event properties. 

The same formula and code is used for producing graviton events. In this case A = A' = 1 is 
used independent of the input values and the djj-body phase space factor is changed to, 

16vr 2 0F T(du + \) 27T-W 2 
A{du) = (2vr)^ T( du - l)T(2 du ) ~* S {n) = ^(fT (?) 

The remaining input arguments are used in the following way, 

du = l + l (8) 
Ac/ = M D (9) 

which implies that in total only two constants have to be changed in order to switch between 
the models. 

A truncation option was also implemented in order to test the validity of the effective theory. 
In order to trust the perturbative calculations of the effective theory the s of the process should 
be smaller than the fundamental mass scale of the theory, Kjj or Mjj. This does not have to 
be the case at the LHC where the collision energy is 14 TeV and, for example, the Md scale at 
present does not have to be larger than about 1 TeV. For this reason the truncation functionality 
suppresses the cross section at s > Afj by a factor Afj/s 2 . This truncation also implies that 
the mass spectrum is suppressed at large values. The truncation effect becomes increasingly 
significant with large djj values, since the mass spectrum then becomes peaked towards higher 
values. This is illustrated in figure 1, which shows the graviton mass spectrum for n = 2 with and 
without truncation (left) as well as for n = 1 and 6 with no truncation (right). The truncation in 
the case n = 2 is using Mo = 2 TeV. In a similar way the truncation effect becomes increasingly 
important with an increased transverse energy requirement of the event selection in the analysis. 

The truncation option was primarily used to verify that the difference between the cross 
section obtained with and without truncation is negligible. This implies that the uncontrolled 
s-region does not contribute significantly. However, in the case of a discrepancy, the trun- 
cated cross section can still be used as a conservative estimate from the region which is under 
perturbative control. 

4 Validation Cross Checks 

In order to validate the different parts of the implementation, the results were checked against 
several similar processes: 

• e+e" -► 7 + G in [2, 3]; 

• pp — »■ 7 + G in [2] ; 

• e+e~ -► Z + G in [8, 13]; 

• e + e~ -> 7 + U in [11]. 
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Figure 1: The graviton mass spectrum for n = 2 with and without truncation (left) as well as for n = 1 
and 6 with no truncation (right). 



Together these should provide sufficient cross checks to validate all parts of the implemented 
process. Most parts are covered by reproducing the results in [2, 3] from the first two processes. 
The third processes verifies the Z specific parts of the spin-2 matrix element. The fourth process 
verifies that the unparticle phase space factors and spin-1 matrix element are correct. The only 
part which is not directly tested is the matrix element in the case of a spin-2 unparticle with 
A A'. On the other hand, when setting A = A' a specific cancellation is required in order to 
reproduce the Z + G matrix element and therefore this is considered to be a good general test. 

The same code was also used to produce the photon results. The photon limit of the Z+G/U 
process was obtained by making the following changes, 



cos 9\y 



m z ^0 (10) 

Q 2 (ii) 

(12) 



9v + 9 a _^ q2 
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where the following coupling conventions were used, g v = — \ + 2 sin 2 9w and g a = —\. The fact 
that the photon limit of the Z matrix element is checked is also a valuable cross check of the 
implement ation . 



5 Z + G/U Production at the LHC 

The potential of using Z + G/U events to confirm a signal observed in the main channels was 
also studied. These events would give rise to di-leptons and large missing transverse energy. 
The most favorable event topology from an experimental point of view is when the leptons are 
either an electron or muon pair. The SM process with the most similar experimental signature 
is ZZ production, where one Z decays into electrons or muons and the other into neutrinos, 
ZZ — > ttvv. For this reason it is interesting to investigate the possible Z + G/U cross section 
relative to this SM ZZ process. 
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The ATLAS ZZ — > ttvv analysis presented in [15] will be used as a reference through this 
study. This analysis predicts an overall inclusive ZZ cross section of ozz = 14.8 pb, which was 
obtained by the MC@NL0 generator and using the CTEQ6M parton distribution functions (PDFs). 
The corresponding cross section obtained by Pythia 8 was 11.4 pb using the CTEQ5L PDFs. 
This appears consistent for a comparison between leading order (LO) and next-to-leading order 
(NLO) results. 

The ZZ events generated by Pythia were then passed through a similar event selection to 
the one used in [15]. This corresponds to the following cuts: 

1) Opposite charged electrons or muons; 

2) p l T > 20 GeV; 

3) \rf\ < 2.5; 

4) \M U ~ 91.2| < 10 GeV; 

5) p T {Z) > 100 GeV. 

The first three requirements ensure that there are two high-p^ electrons or muons inside the 
central part of the experiment. The remaining two cuts imply that the leptons are consistent 
with the decay of an on-shell Z with high pr- After applying this selection to the generator 
level ZZ events from Pythia, a selected cross section of 14 fb was obtained. This value also 
appears reasonable compared to the 10.2 ± 0.2 fb obtained in the full ATLAS analysis, when 
considering that they are calculated at LO and NLO respectively, and that the ATLAS analysis 
includes the full ATLAS detector simulation and reconstruction chain. 

Since the existing limits on Md exclude values between approximately 0.5 and 1.5 TeV the 
following parameter values were used as benchmark points for the Z + G study: 

• n = 1, 2, 3, 6; 

• M D = 2, 2.5, 3 TeV. 

Events generated with, n = 2, = 2 TeV and using the truncation option were passed through 
the ZZ selection described above and a selected cross section of 2.6 fb was obtained. This shows 
that for this particular point the Z + G signal would amount to about 20% of the expected SM 
ZZ — > Itvv signal. 

Figure 2 shows the selected cross section for the different benchmark points and the results 
are shown with and without truncation. One can clearly see that the cross section increases 
with a larger number of extra dimensions, but at the same time an increasing fraction moves 
into the non-perturbative s region. For this reason the effective theory is not trustworthy at 
large n. At small n, on the other hand, the effective theory is fairly valid and the Z + G channel 
could be used to cross check a signal observed in the more conventional channels. This will, 
however, require relatively large amounts of data. Assuming for example 05 = 20 fb of selected 
background 3 in total and a signal of a s = 2.6 fb, a data sample in the order of, 

C = ~ 75 fb' 1 (13) 

3 This background value is assumed since the ATLAS study suggests about 50% more selected ZZ candidates 
than real ZZ events. 
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Figure 2: Selected Z + G cross section as a function of Mo- The results are shown with (dashed line) 
and without (solid line) truncation. 



would be necessary for a 5a discovery. 

The event selection could be improved to increase the sensitivity of the graviton signal. The 
left plot of figure 3 shows the pr{Z) distributions of the SM ZZ — > llvv process and the Z + G 
process with n = 2 and Mjj = 2 TeV. Here a different shape can be seen and this could, 
for example, be used to optimize the selection. It was also observed that the pr(Z) spectrum 
becomes harder with an increasing value of n. A detailed study of the selection should, however, 
include proper detector simulation and event reconstruction. This will not be covered here, but 
would be a natural continuation of this work. 

The same selection was used to investigate the case of spin-1 unparticle emission. Here the 
Z + U signal was generated using the parameter points d\j = 1.4, 1.6 and 2.0 in accordance with 
[11] and A = 1 was kept fixed. Despite the fact that the unparticle mass spectrum shown in 
equation (6) is the same as in the spin-2 case, the spin-1 matrix element implies a much softer 
mjj as well as pt{Z) distribution. The pr(Z) distribution is illustrated in the right plot of figure 
3, which shows spin-1 unparticles with djj = 2 compared to the SM ZZ distribution. For this 
reason the detection efficiency is significantly lower than in the spin-2 case and, even with large 
amounts of data, the prospects seem poor for reaching significantly beyond the LEPII limits [11] 
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Figure 3: Comparison between the pr(Z) spectra from the SM ZZ — > ttvv and the Z + G/U processes 
including truncation. The left plot corresponds to a Z + G signal with n = 2 and Md = 2 TeV. The right 
plot corresponds to a spin-1 Z + U signal with djj =2 and Kjj — 2 TeV. All distributions are based on 
events with two electrons or muons that pass the cuts 1 to 3 and have the same arbitrary normalization. 



using the Z + U channel at the LHC. 

The corresponding results for the spin-2 unparticle scenario with A = A' can be obtained 
simply by scaling the graviton results according to, 



otj A(du) 



A 2 (14) 



oq S'(n) 



For this reason, the study of spin-2 unparticles was focused on deviating effects when A / A'. 
Figure 4 shows the my distributions of events generated with djj = 2 and r[= A'/A] = 1.01 as 
well as r = 1. In the left plot it can be seen that even for this small deviation from A = A', 
a large spike develops in the first bin. This spike is resolved in the right plot and indicates 
an infra-red (IR) divergence of the cross section in the case of A / A'. This is explained to 
some extent by the fact that the terms in the matrix element that cancel when A = A' are all 
proportional to powers of l/mfj. For this reason also the total inclusive cross section increases 
rapidly when the couplings are not equal. Unlike the spin-1 case, the events in the soft mjj tail 
have a significantly harder pr(Z) spectrum than the SM ZZ process. This indicates that the 
Z + U channel could also be a useful complement when investigating scenarios with A / A'. 



6 Conclusions 

Models with extra dimensions have attracted much interest as candidates for potential physics 
beyond the SM. One of the most popular phenomena is real graviton emission and it has been 
pointed out earlier that graviton emission can be seen as a special case of the so-called unparticle 
emission process. This allows the same MC implementation to cover both scenarios and the 
parton level process // — > Z/^f + G/U has been implemented in Pythia 8. In addition, Z + G/U 
production at the LHC has been studied as a complement to the more conventional channels. 
This process was studied at generator level relative to the similar SM process ZZ — > ttvv. In 
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Figure 4: The (spin-2) unparticle mass spectra for djj = 2 with r = 1 and r = 1.01. The histograms are 
arbitrarily normalized. The distributions in the left plot, however, have the correct ratio with respect to 
the total cross sections. The cut at 1 GeV in the right plot is a free parameter of the MC implementation. 

order to get a handle on the experimental effects, a detailed analysis of the ATLAS ZZ — > ttvv 
cross section measurement has been used as a reference. The study indicates that Z + G 
production could be used for small values of n to confirm a signal observed in the main channels. 
It will, however, require a relatively large amount of data. Spin-1 unparticle production at the 
LHC, on the other hand, seems more difficult to observe using the Z + U channel. Spin-2 
unparticle events have the same characteristics as the graviton events if A = A'. When these 
couplings are different an IR divergence is observed in the mu spectrum. This increases the 
total inclusive cross section and the Z + U channel could provide a useful complement also in 
this scenario. 
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Appendix - Implementation in Pythia 8 



The processes ff—*Z + G/U and // — > 7 + G/U are implemented in two semi-internal classes 
called Sigma2f f bar2UZ and Sigma2f f bar2Ug. The classes include a number of standard Pythia 
methods which are called at different points during the event generation. For this reason they 
can be used by the main Pythia 8 library in the same way as for example the Sigma2bg2Hb 
process in the main25.cc example of Pythia version 8.108. 

The constructors takes 7 arguments in order to instantiate an object, 

• int Spin, unparticle spin (1 or 2); 

• bool Trunc, set true in order to truncate the contributing s values; 

• bool Graviton, set true to use graviton specific settings according to eqs. 7, 8 and 9; 

• double du, scale dimension parameter; 

• double Ajj, unparticle renormalization scale; 

• double A, unparticle coupling to SM; 

• double Ratio, the ratio A'/A of the spin-2 matrix element [11]. 

If Graviton = true, the A and Ratio values are overridden and set equal to one. 

The graviton particle code 39 was used both in the case of graviton and unparticle emission. 
The graviton available in Pythia has a Breit-Wigner distributed mass and the underlying Breit- 
Wigner shape can be adjusted in the main program by the following parameters, 

pythia. readString(" 39 :m0 = 50."); 
pythia. readString(" 39 :mWidth = 500."); 
pythia. readString(" 39 :mMin = 1."); 
pythia. readString(" 39 :mMax = 13910."); 

The Breit-Wigner was re-weighted to a flat distribution by using the runBW3 weight in the 
sigmaHatO method. By doing this the differential unparticle cross section in equation 6 en- 
sures a proper mass distribution of the generated final events. 

Program Summary 

Code Location: sask.home.cern.ch/sask/gravunplOO.tgz 
Tested with: PYTHIA version 8.108 

The . tar file contains: 
README 
main-gravunp. cc 
Sigma2ffbar2UZ.h 
Sigma2ffbar2Ug.h 
Sigma2ffbar2 UX. cc 
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